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Abstract 



Thermal Hilbert moment QCD sum rules are used to obtain the temperature dependence of the 
hadronic parameters of charmonium in the vector channel, i.e. the J/ip resonance mass, cou- 
pling (leptonic decay constant), total width, and continuum threshold. The continuum threshold 
sq, which signals the end of the resonance region and the onset of perturbative QCD (PQCD), 
behaves as in all other hadronic channels, i.e. it decreases with increasing temperature until it 
reaches the PQCD threshold sq = 4mQ, with mg the charm quark mass, at T ~ 1.22 Tc. The 
rest of the hadronic parameters behave very differently from those of light-light and heavy-light 
quark systems. The J/ip mass is essentially constant in a wide range of temperatures, while the 
total width grows with temperature up to T ~ 1.04 Tc beyond which it decreases sharply with 
increasing T. The resonance coupling is also initially constant and then begins to increase mono- 
tonically around T ~ Tc. This behaviour of the total width and of the leptonic decay constant 
provides a strong indication that the J/tp resonance might survive beyond the critical tempera- 
ture for deconfinement. 
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1 Introduction 



A successful quantum field theory framework to extract hadronic information from QCD analyt- 
ically is that of QCD sum rules [1]. This technique is based on the Operator Product Expansion 
(OPE) of current correlators beyond perturbation theory, and on Cauchy's theorem in the com- 
plex energy plane (quark-hadron duality) . This program was first extended to finite temperature 
in [2]. It is based on two basic assumptions, (a) that the OPE continues to be valid, with the 
vacuum condensates developing a temperature dependence, and (b) that no thermal singularities 
appear in the complex energy plane, other than on the real axis, i.e. the notion of quark-hadron 
duality also remains valid. Field theory evidence in support of these assumptions was provided 
later in [3]. Numerous applications of QCD sum rules at finite temperature have been made over 
the years [5"[£]) leading to the following scenario for light-light and heavy-light quark hadrons. 
(i) As the temperature increases, hadronically stable particles develop a non-zero width, and 
resonances become broader, diverging at a critical temperature interpreted as the deconfinement 
temperature (Tc). This width is a result of particle absorption in the thermal bath, and resonance 
broadening at finite temperature was first proposed in [7J. (ii) Above the resonance region the 
continuum threshold in hadronic spectral functions, i.e. the onset of perturbative QCD (PQCD), 
decreases monotonically with increasing temperature. In other words, as T — > Tc hadrons melt 
disappearing from the spectrum, which then becomes smooth, (iii) Additional support for this 
picture is provided by the behaviour of hadronic couplings, or leptonic decay constants, which 
approach zero as T — > Tc. Also, hadronic and electromagnetic mean-squared radii diverge at Tc 
[8] indicating deconfinement. On a separate issue, QCD sum rules in the axial-vector channel have 
provided (analytical) evidence for the (almost) equality of the critical temperatures for decon- 
finement and chiral-symmetry restoration [9j . Contrary to this revealing behaviour of widths and 
couplings, the mass does not appear to offer any relevant information about deconfinement. In 
fact, in most cases it increases or decreases, sometimes slightly, with increasing T. Conceptually, 
given either the emergence or the broadening of an existing width, together with its divergence 
at Tc, the concept of mass looses most of its meaning. 

At this stage it must be pointed out that in the framework of QCD sum rules the critical tem- 
perature for deconfinement, referred to above, is only a phenomenological parameter. It is the 
temperature at which the resonance couplings and the continuum threshold approach zero, and 
the widths increase sharply, for light-light and heavy-light quark correlators. Hence, it need not 
coincide numerically with e.g. the critical temperature obtained in lattice QCD [10], which is 
defined differently. In fact, results from QCD sum rules lead to values of Tc somewhat lower than 
those from lattice QCD. Hence, comparisons between different frameworks should be made in 
terms of the dependence of parameters on the ratio T/Tc. 

Turning to heavy-heavy quark hadrons at finite T, in principle one would expect them to be- 
have differently from the light-light and heavy-light quark systems for the following reasons. 
The T dependence of the latter in the PQCD sector is dominated by the spectral function in 
the time-like region, the so-called annihilation term, which is anyway relatively unimportant in 
relation to the light quark condensate contribution. The PQCD spectral function in the space- 
like region (scattering term) is highly suppressed. For heavy-heavy quark systems this is not 
the case; the scattering term becomes increasingly important with increasing temperature while 
the annihilation term only contributes near threshold. In the non-perturbative QCD sector of 
light-light and heavy- light quark correlators, the driving term in the OPE is the light quark 
condensate. This term is responsible for the behaviour of the continuum threshold, as it follows 
that so(T)/so(0) ~ {{qq)) / {qq) [5]-[H1. This scaling relation is rather important as the thermal 
light quark condensate is the order parameter for chiral symmetry restoration. In contrast, for 
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heavy-heavy quark correlators the leading power correction in the OPE is that of the gluon con- 
densate, which has a very different temperature behaviour. This expectation about a different 
temperature behaviour of heavy-heavy quark systems is shown here to hold for the case of the 
J/il). Using Hilbert moment QCD sum rules we obtain the J/tp hadronic parameters at finite 
temperature. With the exception of the continuum threshold, we find a very different behaviour 
from that of light-light and heavy- light quark systems. The continuum threshold, so(T), does 
decrease with increasing T, being driven by the gluon condensate and the PQCD spectral func- 
tion in the space-like region, until it reaches the PQCD threshold sq = where mq is the 
charm quark mass, at T ~ 1.22 Tc. The J /if: mass remains basically constant. The width and 
the coupling are almost independent of T up to T ~ 0.8 Tc where the width begins to increase 
substantially, but then above T ~ 1.04 Tc it starts to decrease sharply, and the coupling increases 
also sharply. This behaviour, which can mostly be traced to that of the PQCD spectral function 
in the space-like region, points to the survival of the J/il) resonance above the deconfinement 
temperature. However, the QCD sum rules have no longer solutions for the hadronic parameters 
once the continuum threshold reaches the value so|mm = -^j/^! there is no longer any support 
for the integrals. Hence, the temperature range explored with this technique does not extend be- 
yond T ~ 1.22 Tc. Given the importance of the PQCD spectral function in the space-like region, 
it should be realized that non-relativistic approaches to charmonium at finite temperature will 
most probably miss this contribution. In fact, the complex energy plane in the non-relativistic 
case would only have one cut along the positive real axis, which would correspond to the time- like 
(annihilation) region of PQCD. The space-like contribution (g^ = {uJ^ — |q|^) < 0) in the form of 
a cut in the energy plane centered at the origin for — |q| < w < |q|, would not be present in the 
non-relativistic case. 



2 Hilbert moment QCD sum rules 



We consider the correlator of the heavy-heavy quark vector current at finite temperature 
Yl^,{q\T) = i j d^xe"^"^ 9{xo) « \[V^ix) , m\ » 

= -{oi^uq^ - qpiqu)^{q^,T) , (l) 

where V^(x) =: Q{x)'y^Q{x) :, and Q{x) is the heavy (charm) quark field. The vacuum to vacuum 
matrix element above is the Gibbs average 

« A-B »= exp{-En/T) < n\A ■ B\n > /Tr{exp{-H/T)) , (2) 

n 

where |n > is any complete set of eigenstates of the (QCD) Hamiltonian. We shall adopt the 
quark-gluon basis, as this allows for the standard QCD sum rule program to be smoothly extended 
to finite temperature [3]. 

The imaginary part of the vector correlator in perturbative QCD (PQCD) at finite temperature, 
Im n(g^,T), involves two pieces, one in the time-like region {q^ > 4m^), Im na{q^,T), which 
survives at T=0, and one in the space-like region (q^ <0), Im Us{q'^,T), which vanishes at T=0. 
A straightforward calculation in the time-like region, to leading order in PQCD, gives 
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where = l — Arnq/tf', uiq is the heavy quark mass, q'^ = co'^ — c^ > 4mQ, and npiz) = (l+e^)~^ 
is the Fermi thermal function. In the rest frame of the thermal bath, |q| — > 0, the above result 
reduces to 

- Im Uaiio, T) = v{3 - v^) [1 - 2nF{uj/2T)] e{oj - 2mQ) . (4) 

The quark mass is assumed independent of T, which is a good approximation for temperatures 
below 200 MeV [llj . As is customary in all QCD sum rule analyses at finite temperature, only the 
leading order in the strong coupling will be considered here. One reason is that the temperature 
introduces an additional scale, and this problem is not yet fully understood. More importantly, 
though, results in this framework are not intended to be of high precission, as the T-dependence 
of hadronic parameters will probably never be measured with great accuracy (some, like e.g. 
the leptonic decay couplings, may not even be measured at all at finite T). For this reason one 
normally determines in this framework the ratio of hadronic parameters at finite and at zero T 
as a function of T/Tc- 



The calculation of the PQCD piece in the space-like region, the so-called scattering term, is more 
involved as the limit |q| ^ must be taken with extreme care. In fact, in the complex energy 
plane, and in the space-like region the correlator n(g-^), Eq.(l), has a cut centered at the origin 
and extending between oj = — |q| and u; = |q|. In the rest frame of the thermal bath this cut 
shrinks to zero and produces a delta function in the imaginary part of n(g^). Details of 

this calculation are left for the Appendix; the result is 
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To complete the evaluation of the vector correlator in QCD, we will add later the leading power 
correction in the OPE, which in this case is given in terms of the gluon condensate ((Ola^C^lO)). 
Turning to the hadronic representation of the vector correlator we shall, as usual, parametrize it 
in terms of the ground state resonance, i.e. the J /tp-, followed by a continuum given by PQCD 
after a threshold sq > My. This ansatz is even a better approximation at finite temperature, as 
sq(T) decreases monotonically with increasing T. Considering first the zero width approximation, 
the hadronic spectral function is given by 

- ImU{s,T)\HAD = -ImU{s,T)\iiES G{so - s) + -ImU{s,T)\pQCD 0{s - So) 

TT TT TT 

= 2f^{T)6{s-M^(T)) +-ImU{s,T)ae{s-so) , (6) 

TT 

where s = = to'^ — , and the leptonic decay constant is defined as 

< 0\V^{0)\V{k) >= V2 Mv fv • (7) 
Next, considering a finite (total) width the following replacement will be understood 

Sis - mUt)) =^ const _ Ml^iTniT) ' 

where the constant is fixed by requiring equality of areas, e.g. if the integration is in the interval 
(0 — oo) then const = My{T)Ty{T)/TT. To complete the hadronic parametrization one needs 
to consider the hadronic scattering term due to the current scattering off heavy-light quark 
pseudoscalar mesons (D-mesons). This contribution is given by [2] 
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where nsiz) = (1 — e^)^^ is the Bose thermal function. As shown in [6] heavy-Hght pseudoscalar 
mesons deconfine at the critical temperature, which in that application was Tc — 100 — 110 MeV. 
Since, a posteriori, the critical temperature for the J/'0 is much higher, this scattering term should 
not contribute. More importantly, though, it is easy to see by comparing Eq.(9) with Eq.(5) that 
this hadronic scattering term is exponentially suppressed; it is, in fact, two to three orders of 
magnitude smaller than the QCD counterpart in the wide range of temperatures explored here, 
to wit. The ratio R of this hadronic scattering term, Eq. (9), and its QCD counterpart, Eq.(5), 
is i? ~ 10-3 at r = 100 MeV, and R ~ 10"^ at T = 160 MeV. 



The correlation function n(g^,T), Eq.(l), satisfies a once subtracted dispersion relation. To 
eliminate the subtraction one can use Hilbert moments, i.e. 



where N = 1,2,..., and Q'^ > is an external four-momentum squared, to be considered as a 
free parameter as discussed below. Using Cauchy's theorem in the complex s-plane, which is 
equivalent to invoking quark-hadron duality, the Hilbert moments become Finite Energy QCD 
sum rules (FESR), i.e. 

m{Q\T)\nEs = m{Q^T)\QCD , (u) 

where 

1 /"^o^^) ds 

'Pn{Q\T)\res ^ - ——-^-^ImIi{s,T)\RES, (12) 



1 f^o(T) 



1 ds 

+ ^ (g + Q2)jV+i IrnUs{s,T)+ipN{Q^T)\NP , (13) 

and ImIl{s,T)\fiES is given by the first term in Eq.(6) modified in finite-width according to 
Eq.(8), and the PQCD spectral functions are given by Eqs.(4) and (5). The non-perturbative 
term corresponding to the dimension d=4 in the OPE is given by [T] 

^-(Q'^n-^ = -i(4^(rT^^(^+''-i^+^'^, 

2^iV(jV + l)^(iV + 2)(iV + 3)(iV-l)! 

(2Af + 5)(2Af + 3)!! ^ ' 

where -P(a, 6, c, z) is the hypergeometric function, , ^ = p = and 
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The gluon condensate ((^C^)) at low temperatures has been calculated in chiral perturbation 
theory [12] with the following result 
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where Np is the number of quark flavours, and Ap 
this can be written as 
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400 MeV. To a good approximation 
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Due to this thermal behaviour the gluon condensate, in this framework, remains essentially con- 
stant up to temperatures close to Tc — 100 MeV, after which it decreases very sharply. In order 
to go beyond the low temperature regime of chiral perturbation theory, lattice QCD provides the 
right tool. A good approximation to results in this framework pLSj is given by the expression 
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Figure 2: The ratio My(r)/My(0) as a function of T/Tc- This ratio is basically the same in 
zero-width as in finite width. 



where T* ~ 150 MeV is the breakpoint temperature where the condensate begins to decrease 
appreciably, and 250 MeV is the temperature at which {{^G'^))rp^ = 0. 

Returning to the dependence of the Hilbert moments, Eq.(lO), it has been customary in many 
analyses of charmonium at T=0 to take = 0, but the case > has also been advocated 
|14j . At finite temperature, though, due to the singular behaviour of the space-like PQCD term, 
Eq.(5), one is compelled to take > 0. Since all hadronic parameters will be normalized to 
their values at T=0, and we are only interested here in their temperature behaviour, we shall 
fix as well as so{0) from the experimental values of the mass, the coupling, and the width 
at T=0. At finite temperature there are additional contributions to the OPE in the form of 
non-diagonal (Lorentz non-invariant) condensates. In the case of non-gluonic operators they are 
highly suppressed [6], [15] so that they can be safely ignored. A gluonic twist-two term in the 
OPE has been considered in [16], and computed on the lattice in [TTj. Using this information we 
find that the non-perturbative QCD moments, ifj^{Q'^,T)\j\[p, change as follows 



where 



1 + i 
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with Go(r) ^ G2 »^ (^^(0) = (0.05 ± 0.02) GeV^ p^, G2(T) ~ -lO^^ GeV^ [17], in the 

range of temperatures considered here, and the parameters ^ and p = (1/2 — 1) have been defined 
after Eq.(14). For small values of N (N = 1 — 3) the second term in brackets above is at the level 
of a couple of percent, while for larger values of A'^ it becomes negligible. The correction term 
is small and essentially independent of Q^, i.e. 5n{Q'^,T) ~ (2 — 6)% in the temperature range 
considered here, and will play no appreciable role in the results, as will be discussed later. 



3 Results 



We begin by determining sq and at T=0 from the moments, Eq.(ll), and using as input 
the experimental values [l9] My = 3.097 GeV, fv = 196 MeV, and Ty = 93.2 keV, as weh as 
niQ = 1.3 GeV, and [18] (0|^G^|0) ~ 5 x 10~'^GeV^. In the zero-width approximation one finds 
from Eq.(12) that 

^i{Q'^)\res _ f2{Q^)\REs .^n 



4000 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 

T/Tc 

Figure 3: The ratio ry(T)/ry(0) as a function of T/T^. 
Given the extremely small total width of the J/^[^ it turns out that the above relation also holds 



with extreme accuracy in finite width. Using Eq.(ll) this leads to 



^i{Q'^)\qcd ^ ^2{Q'^)\qcd .22^ 

^2{Q^)\qCD <^3(Q2)|qcd' ^ ^ 

which depends only on the two unknowns so and Q^, and provides the first equation to determine 
this pair of parameters. The second equation can be e.g. Eq.(ll) with = 1. In this way 
we find that sq = 11.64 GeV^, and = 10 GeV^ reproduce the experimental values of the 
mass, coupling, and width of J/ij) within less than 1%. We have checked that equally reasonable 
agreement for the hadronic parameters at T = can be achieved using larger values of Q^, up to 
~ 20 GeV^, as well as larger values of N , up to say ~ 10. These different choices do not 
produce any qualitative change at finite T, other than some 10 % shift of the critical temperature, 
upwards for larger Q^, and downwards for larger N. This whole set of hadronic parameters will 
then be used to normalize the corresponding parameters at finite temperature. The latter are 
obtained as follows. 

The continuum threshold sq{T) is obtained again using Eq.(22), as Eq.(21) continues to hold to 




Figure 4: The ratio /y(T)//v(0) as a function of T/T^. 



great accuracy in finite width, even if the width were to increase with temperature by three orders 
of magnitude, say from ry(0) ~ 0.1 MeV to Vv{T) ~ 300 MeV. The result for so(T) normalized 
to so(0) is shown in Fig.l as a function of T/Tc, where Tc ^ 140 MeV is the temperature where 
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so(T) begins to deviate appreciably from so(0). Above this temperature so(T) decreases sharply 
until it approaches the threshold Sth = 4mQ at a temperature T ~ 1.22 Tc, with this threshold 
being the minimum possible value for so{T). It must be stressed that the temperature dependence 
of the continuum threshold is rather important in order to avoid inconsistencies, and therefore it 
should not be ignored. In the case of light-light and heavy-light quark systems so{T) is related 
to the thermal light quark condensate. This is not the case for heavy-heavy quark hadrons, 
so that sq{T) must be determined from the Hilbert moments together with the other hadronic 
parameters. Next, the mass, in the zero-width approximation, follows from the ratio 

MQ^t)\res ^ MQ^T)\QCD 

MQ^,T)\res MQ^T)\qcd 

The result for My(r)/My(0) as a function of T/Tc is shown in Fig.2. The mass remains basically 
constant until T Tc where it decreases slightly by a few percent. Beyond T ~ 1.1 Tc the whole 
program of determining the temperature behaviour of the hadronic parameters breaks down. 
The simple reason for this is that at this temperature so(T) has left the resonance behind, i.e. 
Sth < My. Hence, there is no longer a support for the Hilbert moments in the hadronic sector. 
We have verified that this result for Mv'(T)/My(0) remains essentially the same in finite width. 
Hence, in determining the coupling and the width we shall keep the mass independent of the 
temperature. Next, the width can be determined from the ratio 

>fl{Q\T)\RES _ ipi{Q\T)\QCD 

The result for Tv{T) /Ty{{)) as a function of T/Tc is shown in Fig. 3; the peak is reached at 
T ~ Tc- Finally, the leptonic decay constant can be determined from any single moment, e.g. 
from 

MQ\T)\res = MQ'',T)\qcd. (25) 

The ratio fv{T)/ fv{0) as a function of T/Tc is shown in Fig. 4. The rise of the coupling be- 
yond Tc is the result of a detailed balance between different contributions. With increasing T 
the continuum threshold, sq{T), decreases and so does the PQCD annihilation moment, while 
the QCD scattering moment increases approaching the annihilation moment at T ~ 160 MeV; 
all the while, the non perturbative moment remains approximately constant and negative, and 
numerically comparable to the other two terms at that temperature, thus mostly canceling the an- 
nihilation moment there. The hadronic moment dividing the QCD moment is mostly a decreasing 
function of the temperature due to the decrease of so (T) and to the increase of the width, leading 
finally to an increase in the coupling. This behaviour of the coupling and width of the J/ remains 
qualitatively stable against changes in N and in the moments, for ~ 10 — 20 GeV^, and 
~ 1 — 10. Numerically, higher values of tend to increase slightly the critical temperature, 
and higher values of tend to reduce it. We have explicitly verified that the hadronic scattering 
moment can be safely ignored. Both the width and the coupling can only be determined up to 
Tf ~ 1.1 Tc beyond which sq(T) < My{T). The temperature behaviour of the width and the 
coupling thus obtained strongly suggests the survival of the J/ip above the critical temperature 
for deconfinement. 
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(23) 



(24) 



4 Discussion and conclusions 



In this paper we have determined the temperature behaviour of the hadronic parameters of the 
J/ijj resonance using Hilbert moment QCD sum rules. The hadronic parameters are not only the 
mass, coupling (leptonic decay constant), and total width, but also the continuum threshold sq 
which signals the onset of PQCD. This threshold at finite temperature was first introduced in [2], 
where it was interpreted as a phenomenological parameter signaling deconfinement. Later, it was 
related to chiral symmetry restoration [9j, and a scaling law was found associating sq{T) to the 
thermal (light) quark condensate [8]. Hence, this parameter constitutes an essential part of any 
QCD sum rule analysis at finite T, and ignoring the temperature dependence of sq would lead to 
inconsistencies. The results of our analysis at low temperatures show that the J/ijj behaves as 
other light-light and heavy-light quark resonances, i.e. the continuum threshold and the coupling 
decrease and the width increases with increasing T. However, as T approaches Tc this behaviour 
changes so that the coupling increases and the width decreases. The PQCD spectral function 
in the space-like region (scattering term), as well as the gluon condensate, are responsible for 
this new scenario which strongly suggests the survival of the J/ij: beyond deconfinement. This 
is a feature unique to heavy- heavy quark systems, as this space- like contribution is negligible in 
light-light and even heavy-light quark channels [6j. While this PQCD term increases monoton- 
ically in importance with increasing T, the QCD sum rule framework necessarily breaks down 
at the temperature where sq{T) reaches the minimum value so| mm = 

M'^jj^. In fact, since the 

hadronic mass is essentially independent of T, once sq{T) < Mj^^ there is no further support 
for the Hilbert integrals. In spite of this, one could continue to explore the hadronic parameter 
space by reducing arbitrarily the value of Mj/^ at T = 0. In this case, one would find that 
the behaviour of the coupling and width persists at higher temperatures. Our results are not 
necessarily in conflict with a recent QCD sum rule analysis [T7], where charmonium was found 
to survive beyond T^. The reason is that in [T7j the scattering term was not taken into account, 
and the temperature dependence of the continuum threshold was not considered. Both these 
features are essential for a consistent QCD sum rule analysis. Studies of charmonium at finite T 
in non-relativistic frameworks [20] -[21] also lead to charmonium melting at some critical tempera- 
ture. However, in the non-relativistic regime the complex energy plane has only a right-hand cut 
extending from zero to infinity, which corresponds to the time-like case in the relativistic domain. 
The space- like cut at finite T centered at the origin in the complex energy plane, and covering 
the range — |q| < to < -|-|q|, is a genuine relativistic effect. The results and conclusions of the 
present analysis, though, are in agreement with those obtained from QCD numerical simulations 
on the lattice |10j . 



5 Appendix 

In this appendix we outline the calculation of the imaginary part of the correlator, Eq.(l), in 
PQCD and in the space-like region, as there are some incorrect results in the literature. To 
one-loop order in PQCD this is given by [2] 




(26) 
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where 



(27) 



Performing the change of variable y = -^^j the Fermi factor Ani? in Eq.(26) becomes the total 
derivative 



iXuF = — — nFyV) 



(28) 



After integration by parts, Eq.(26) becomes 



ImYl,{q^,T) 
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Finally, taking the limit w — > 0, followed by |q| 0, and using the expression 



lim lim ( - — -3- | = - 6 (uj'^) , 
IqHO V|q|V 3 ^ ^ ' 



(30) 



gives the final result for the PQCD spectral function in the space-like region, i.e. 



— ImIls{uj,T) = —2 uiq 5(0;^ 



vr 
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2 T 

/ ynF{y)dy 



(31) 
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